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Q\ ■ Abstract 



Passive scalar transport in a spatially-extended background of roll convection is 
considered in the time-periodic regime. The latter arises due to the even oscillatory 
instability of the cell lateral boundary, here accounted for by sinusoidal oscillations 
of frequency uj. By varying the latter parameter, the strength of anticorrelated re- 
gions of the velocity field can be controled and the conditions under which either an 
enhancement or a reduction of transport takes place can be created. Such two ubiq- 
uitous regimes are triggered by a small-scale(random) velocity field superimposed 
to the recirculating background. The crucial point is played by the dependence of 
Lagrangian trajectories on the statistical properties of the small-scale velocity field, 
e.g. its correlation time or its energy. 

PACS 47.27Qb - Turbulent diffusion. 

Transport in turbulent flows with recirculation is a problem of great interest both in 
the atmosphere and in the ocean [|l], |2|, |3|. In the atmosphere, the so-called horizontal roll 
vortices or, briefly, rolls are the paradigm of recirculating pattern. They can be considered 
as a manifestation of Benard-Rayleigh convection occurring when conditions of combined 
surface heating and strong winds take place in the atmospheric boundary layer. Their 
depth equals the mixing layer and the ratio of lateral to vertical dimensions for a roll 
pair is about 3:1 (see, e.g., Ref. ||). In the following, such type of flow configuration will 
be the background where particle tracers (i.e. a passive scalar) will be plunged and their 
statistics investigated. 



An important characteristic of atmospheric rolls is their large spatial extension which 
makes possible the description of the related dispersion problem through an effective 
diffusion equation, i.e. a Fick equation for the large-scale slow-varying passive scalar con- 
centration where the molecular diffusivity is replaced by an enhanced (eddy-) diffusivity. 
The existence of an asymptotic diffusive regime for the large scale concentration can be 
rigorously proved by using, e.g., multiscale techniques ||. A simple Lagrangian interpre- 
tation of the asymptotic diffusive regimes is based on "central limit" argumentations. For 
Lagrangian chaotic flows, velocity correlations decay fastly and, as a consequence, the 
particle displacement, 5x(t), at the time t, is the results of the sum of almost indepen- 
dent advecting contributions. The result is that Sx(t) undergoes Brownian motion when 
observed on times larger than the typical velocity correlation time. 

In this asymptotic framework, the effect on the dispersion process of the small-scale 
components of the velocity field is the renormalization of the effective diffusion coefficient. 
Notice that, as pointed out in Ref. ||, an eddy-diffusivity based description should not 
be possible in the presence of finite-size domains with a small number of recirculations 
where, necessarily, the asymptotic regime might not be reached and the dynamics is 
governed by transient behaviors 0. Interesting studies on this regimes can be found, e.g., 
in Refs. [g g 0. 

For the system we are going to investigate, the characteristic length of the organized 
array of cells is smaller than the size of the domain. As a consequence, the eddy-diffusivity 
tensor, defined through the following asymptotic limit, 



D% = Umt^ — dx^t) - x a (0)}[xp(t) - ^(0)]>, (1) 

with x(t) being the particle position at the time t, and (•) the average over an ensemble 
of tracer particles, turns out to be a well-defined mathematical quantity. 

Atmospheric rolls show a wide range of regimes ranging from (almost) time indepen- 
dent to turbulent flow. The related dispersion phenomena are clearly strongly influenced 
by these different regimes and, as a consequence, transport rates vary over a wide range. 



An intermediate regime attracting considerable attention both theoretically [11 , experi- 



mentally Jl2[ and numerically |13] is the time-periodic regime, where the transport process 



is dominated by advection of tracer particles across the lateral boundary. Such regime 
will be the main concern of the present Letter. Specifically, the main question addressed 
here concerns the role of the small-scale (not explicitly resolved) components of the ve- 
locity field on the large-scale transport. We shall show that the superposition of a colored 
(random) noise velocity (that can be though as associated to a small-scale turbulent mo- 
tion with nonvanishing memory) to the convective (deterministic) background strongly 
affects the transport process: either an enhanced or a reduced (with respect to the white 
case) eddy-diffusion may occur, depending on the frequency of the lateral roll oscillations. 



The interference mechanism, recently proposed in Refs. |14|, 15] for the simple, idealized 



parallel flow, is identified here as the responsible of this twofold behavior. 



Our two-dimensional model for the roll-convection follows Ref. [TI|. Specifically, the 
convective flow is defined by the following stream function : 

ip(x,y, t) = ip sm[k x (x + Bsmut)] sin(k y y) , (2) 

where y ranges from to L y = 2ir/k y . The stream function (||D describes single-mode, 
two-dimensional convection with rigid boundary condition, where the even oscillatory 
instability is accounted for by the term B sin ut, representing the lateral oscillation of 
the roll. In Ref. pi, a quantitative comparison of the behavior in this flow with the 



experimental data has shown that the basic mechanisms of convective transport are well 
captured by the expression (H). 

The periodicity of the cell along the x-axis is denoted by L (L = 2Tc/k x ) while its depth 
(along the y— direction) is L/3 (i.e. k y = 3k x ). The amplitude, B, of the roll oscillations 
is assumed ~ 0.13L. The dimensionless parameter controlling the dynamics is e = ou/ur, 
ujr = k x k y ipQ being the characteristic frequency of particle oscillations inside the cell. The 
two limiting regimes e <C 1 and e ^> 1 have been investigated analytically in Ref. Jll] 



to obtain expressions of the eddy-diffusivity in the limit of zero molecular diffusivity. 
Here, we shall concentrate on the behavior for a wide range of e and in the presence 
of diffusivity. The investigation is however not accessible by analytical techniques. In 
order to evaluate eddy-diffusivities, we have therefore decided to perform Monte Carlo 
numerical simulations of the Langevin equation 

^l = v(x(t)) + v>(t) . (3) 

The velocity field v(x(t)) is incompressible, and related to the stream-function ([|) through 
the usual relations v x = d y ip, v y = —d x ip. The noise term v'(t) is a Gaussian, zero-mean 
random process with the colored-noise correlation function: 

HW#')> = -V^ , (4) 

T 

where Do can be though as the (isotropic) eddy-diffusivity arising from the smallest (not 
explicitly resolved) scales of turbulent motion and r is their correlation time. Notice that 
the white-in-time correlation function is obtained by taking the limit r —> 0. 

From the Langevin equation (Q) the expression ([[]) for the eddy diffusivity can be 
easily rewritten in terms of the Lagrangian autocorrelation C a p(t) = (v a (x(t))vp(x(0))): 

D^ = D o 6 a0 + dtC af3 (t). (5) 

</ 



The role played by anticorrelated regions of the velocity field (i.e. regions where C a p < 
in Eq. ([5])) on the large-scale transport has been investigated in Refs. |TJ], T5\ for the class 



of parallel flows. Two different regimes of transport may occur depending on the extension 
of such regions. Specifically, when anticorrelated regions are sufficiently extended (we 
denote such regime with the label EAR), an increasing r (for a fixed, small, Dq) causes a 
reduction of transport (with respect to r = 0) while an increasing Dq (for a fixed, small, 
r) leads to an enhancement of transport (with respect to Dq = 0). 

The scenario is opposite for anticorrelated regions weak enough (hereafter regime 
WAR), that is, D leads to transport reduction while r to transport enhancement. 
We briefly recall the basic mechanisms characterizing the above two different regimes. 
The first mechanism works to increase the Lagrangian correlation time (and thus eddy- 
diffusivities) : this is due to the fact that r makes the particles of diffusing substance 
forget their past less rapidly than in the case r = 0. Thus, the autocorrelation function 
in ([Sp decays less rapidly than in the withe-noise case. This implies an increasing weight 
of regions where the velocity is strongly (positively) correlated and, as an immediate 
consequence, an increasing eddy-diffusivity. 

The second mechanism arises for flows with closed streamlines and it is associated to the 
presence of anticorrelated regions of the velocity field. The correlation time, r, is now 
working to increase effects of trapping due to the anticorrelated zones where the velocity 
is weak. This means that regions where the velocity is anticorrelated give an enhanced 
(again with respect to the withe- noise case) contribution to the time- integral (||). The 
contribution of anticorrelated regions to the time-integral being negative, a reduction of 
diffusion occurs. 

Our main aim here is to show that in the presence of roll convection the aforesaid 
two mechanisms are relevant and work in competition thus governing the large-scale 
transport. The frequency, uj, of the lateral roll oscillation is identified here as one of the 
parameters controlling the crossover from the WAR and the EAR regimes. Notice that, 



unlike Refs. [[LJ, [15[], our control parameter, uj, is not trivially related to the extension of 



anticorrelated regions. The relation is intrinsic and selected by the dynamics. Indeed, by 



varying uj, it is possible to synchronize JT6| the frequency (of order of ujr) of particles inside 
the cell with the frequency, uj, of the lateral roll oscillation. Due to this synchronization 
mechanism, the eddy-diffusivity as a function of uj can have maxima (when oscillations are 
in phase) or minima (when oscillations are in phase opposition). Moreover from Eq. (|j) it 
results that maxima (minima) of diffusion are associated to the flow configurations with 
the weakest (strongest) anticorrelated regions. 

In order to evaluate the component of the eddy-diffusivity along the direction (e.g., 
the x-axis) of the lateral roll oscillation as a function of uj, numerical integration of Eq. @ 
has been made by using a second-order Runge-Kutta scheme and then performing a linear 
fit of ([x(t) — a;(0)] 2 ) vs t. Averages are made over different realizations and performed 
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Figure 1: Behavior of the eddy-diffusivity (component along the direction of the lateral 
roll oscillations) vs the lateral roll oscillation frequency. Full line: r = 0; dotted line: 
T~/t R = 0.24; long-dashed line: r/t R = 0.48; dot-dashed line: r/t R = 0.95. 

by uniformly distributing 10 6 particles in the basic periodic cell. The system evolution 
has been computed up to times 10 4 t R , where t R = 2tt/u r . 

The x-component, D E , of the eddy-diffusivity us the frequency, u, of the roll oscillation 
is shown in Fig. 1 for Dq/i^o = 5 x 10~ 3 and different values of r: r = (full line), r/t R = 
0.24 (dotted line), r/t R = 0.48 (long-dashed line) and r/t R = 0.95 (dot-dashed line). A 
few comments are in order. Eddy-diffusivity shows maxima originated from the resonance 
between the lateral roll oscillation frequency and the characteristic frequencies of the 
particle motion. Moreover, the effect of r on the shape of the peaks is twofold : firstly, the 
variation of r causes a shifting of maxima positions and, secondly, the larger (smaller) r, 
the higher (lower) the peaks. The first feature (particularly evident for u> corresponding 
to the highest peak) suggests that the correlation time, r, acts to renormalize the large 
scale velocity. The result is that an enhanced (with respect to the white-in-time case) 
convective velocity governs the transport and, in order to have resonance when increasing 
r, the roll oscillation frequency must thus follow the increasing velocity. The renormalizing 
effect of t has been identified perturbatively in Ref. |T7[ for small r. Our results seem to 
suggest a generalization for finite r. 

Concerning the second effect played by r, this means that peaks and valley of the eddy- 
diffusivity profile are associated to regions of type WAR and EAR, respectively. In terms 
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Figure 2: Behavior of the eddy-diffusivity vs the correlation time of the small-scale velocity 
v'. On the left: u/ur=0.67, corresponding to a local minimum in the eddy-diffusivity 
profile of Fig. 1; on the right: u/ur=1.1, corresponding to the highest peak (see again 
Fig. 1). In both cases, D /i/) = 5 x 10~ 3 . 



of the two mechanisms associated to such regions, the first is the winner for value of 
uj corresponding to the peaks in the eddy-diffusivity profile (where the contribution of 
anticorrelated regions is weak) , while the second dominates for values of u corresponding 
to local minima in the eddy-diffusivity profile (where the weight of anticorrelated regions 
is strong). This can be easily seen also from Fig. 2 where behaviors of the eddy-diffusivity 
as a function of r are shown for two different value of u: oj/uor = 0.67 (on the left) and 
uj/ujr = 1.1 (on the right). The former value corresponds to a local minimum in the 
eddy-diffusivity profile, while the latter to the highest peak (see Fig. 1). 

It is remarkable that the scenario above described is opposite when fixing r and varying 
the (small-scale) bare eddy-diffusivity Dq. This can be easily observed from Fig. 3 where 
behaviors of the eddy-diffusivity are now shown as a function of Dq, keeping r fixed and 
equal to zero. As in Fig. 2, on the left we have lo/uir = 0.67 while, on the right, uj/ujr = 1.1. 
The physical reason of such behavior can be easily grasped from the aforesaid mechanisms 
but now recalling the fact that an increasing D makes the particles of diffusing substance 
forget their past more rapidly (rather than less rapidly as it happens by increasing r). 
This effect causes a reduction in the transport. On the other hand, concerning the second 
mechanism, trapping due to anticorrelated regions is now less effective. Indeed, leaving 
the region of trapping is easier when increasing Dq. This fact leads to a reduction of the 
weight of the negative contribution to the time-integral (|5|) giving the eddy-diffusivity 
and, as a consequence, transport is enhanced. 
The final result is thus a complete symmetry between the following operation: increasing 
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Figure 3: As in Fig. 2 but the eddy-diffusivity is plotted vs D for a fixed value of the 
correlation time: r = 0. 



r (for a fixed Dq) and decreasing Dq (for a fixed r). Large-scale transport seems thus 
controlled by a parameter oc Dq/t. Roughly speaking, this is easily understood if we 
observe that when increasing r the particle motion along the Lagrangian trajectories 
becomes more and more coherent. Conversingly, coherence becomes lost when, for a 
fixed r, we decrease Dq. In this sense, the dependence of Lagrangian trajectories on the 
statistical properties of the small scale velocity is crucial in this problem. 

In conclusion, two ubiquitous different regimes of transport have been identified here 
as relevant in the time-periodic roll circulation. A key role to select such two regimes 
of transport enhancement/reduction is due to the interplay between the even oscillatory 
instability of the cell and the statistical properties of the small-scale (random) velocity 
(e.g. their correlation time or their energy). Specifically, in the model here considered, 
the even instability is accounted for by a sinusoidal lateral boundary oscillation with 
frequency uo, while small-scale velocity activity is described by a bare diffusivity, Dq, 
within a Gaussian, zero-mean random process with correlation time r. When varying uj, 
the eddy diffusivity profile appears very structured with sharp peaks separated by evident 
valley. Peaks turn out to be associated to transport enhancement when increasing r (for 
a fixed Dq) or, conversingly, when reducing Dq (for a fixed r). The situation is opposite 
for values of u corresponding to minima in the eddy-diffusivity profile. The physical key 
role is played by synchronization mechanisms (from which the structured eddy-diffusivity 
profile arises) and by the strength of the anticorrelated regions of the velocity field, the 
weight of which in the time-integral giving the eddy-diffusivity can be either enhanced 
(by reducing Dq/t) or reduced (by increasing Dq/t), thus affecting in different ways the 
large-scale transport. 
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